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Abstract 

Heavy colored scalar particles appear in a variety of new physics (NP) models and could be 
produced at the Large Hadron Collider (LHC). Knowing the total production cross section is im- 
portant for searching for these states and establishing bounds on their masses and couplings. Using 
soft-collinear effective theory, we derive a factorization theorem for the process pp — > SX, where 
S is a color-octet scalar, that is applicable to any NP model provided the dominant production 
mechanism is gluon-gluon fusion. The factorized result for the inclusive cross section is similar to 
that for the Standard Model Higgs production, however, differences arise due to color exchange 
between initial and final states. We provide formulae for the total cross section with large (par- 
tonic) threshold logarithms resummed to next-to-leading logarithm (NLL) accuracy. The resulting 
-KT-factors are similar to those found in Higgs production. We apply our formalism to the Manohar- 
Wise model and find that the NLL cross section is roughly 2 times (3 times) as large as the leading 
order cross section for a color-octet scalar of mass of 500 GeV (3 TeV). A similar enhancement 
should appear in any NP model with color-octet scalars. 
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Discovering the Higgs particle and the mechanism of electroweak symmetry breaking 
is one of the major goals of the Large Hadron Collider (LHC). It is well-known that the 
main Higgs production mechanism is the gluon-gluon fusion process. An important issue in 
determining the total production cross section is the large perturbative corrections in the 
threshold region, defined by z — > 1, where z = m 2 H /s, where rnn is the Higgs mass and s is 
the partonic center of mass energy squared. The leading corrections are enhanced by factors 
of log(l — z)/(l — z) and invalidate fixed order perturbation theory in the threshold region. 
These corrections can significantly affect the normalization of the total cross section even 
though the total cross section receives contributions from a range of s Q, 3BB- In the 
threshold region, the inclusive scattering cross section a(pp — > HX) can be factorized (at 
leading twist) into a hard part, soft part, and parton distribution function (PDF) of gluons 
inside the proton 0], and the renormalization group equations (RGE) for these parts can 
be used to resum the large threhold corrections. For the Higgs production cross section the 
calculations up to the next-to-next-leading logarithm (NNLL) accuracy have already been 
performed (5), [(J and give total cross section about three times bigger than predicted at 
leading order. 

Obviously, properly incorporating these effects will be important for other heavy particles 
predicted in theories of New Physics (NP) that may be observed at the LHC. In this paper 
we will focus on the production of a heavy color-octet scalar, which appears in a number of 
NP models such as grand unified theories , supersymmetric theories |10l. 111!]. Pati- 

Salam unification 12j, |l3j, chiral colo r hj| , and topcolor 15]. We will calculate the cross 
section for the Manohar-Wise model [16| of color-octet scalars which is consistent with the 
principle of Minimal Flavor Violation (MFV) 17|,ll8|]. Performing the resummation for color- 
octet scalars is very similar to the resummation for Higgs [l], y, 0]. Similar resummations 
have been performed for squark-antisquark and gluino-pair production cross sections in 



Refs. 19, 20 
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23 to derive 



We will use Soft-Collinear Effective Theory (SCET) |2J, 
a factorization theorem for the cross section. Gluon-gluon fusion cross sections in the full 
theory are matched onto SCET operators. The matching coefficients for these operators 
will differ between various models, but the structure of these operators is universal. The 
cross section computed with the SCET operators factors into correlation functions of SCET 
collinear and soft fields, and renormalization group equations for these correlation functions 
can be used to perform the resummation of the threshold logarithms. This resummation 
procedure is independent of the NP model. 

Before we discuss the factorization theorem, we will describe some details of the color- 
octet scalar model we will be focusing on. The principle of MFV requires that the Yukawa 
couplings of the color-octet scalars be proportional to the Yukawa matrices in the Standard 
Model. The Standard Model Yukawa couplings are 

£ = - gf 3 d Ri H^Q Lj - gVu m eH*Q Lj + h.c. , (1) 

where Ql is the doublet of lefthanded quarks, ur and cLr are the righthanded up and down 
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quarks, respectively, i and j are flavor indices, and H is the Higgs doublet: 
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(2) 

If the color-octet scalar Yukawa couplings are 

C = - \°d Rl S a ^T a Q Ll - \Vu m eS a *T a Q Lj , +h.c. . (3) 

then the MFV hypothesis requires A^' D = r]u,Dgfj C where t)u,d are constants. This elimi- 
nates tree-level flavor changing neutral currents and ensures that experimental constraints 
from flavor physics are not violated. Note that this also implies that the S a couple most 
strongly to the third generation of quarks. The color-octet scalars have gauge couplings to 
gluons and a gauge invariant mass term 

Cqcd = ~S a (D 2 rS c - \m 2 s S a S a , (4) 

where D ac = dj ac + gf abc A*. Finally, there is a scalar potential for the S a and the Higgs 



that can be found in Ref. [16J . 

Though it may seem ad hoc to impose A^' D oc g^' D , this can arise naturally in certain 



models. For example, consider the chiral color model of Ref. [14J]. In this model, the 
gauge group is enlarged to SU(3)l x SU(3)r x SU(2)l x U(1)y, and the chiral color group 
SU(3)l x SU(3)ji breaks down to SU(3) C at some high energy scale. If the righthanded 
quarks are placed in the (1,3,1) representation of SU(3)l x SU(3)r x SU(2) l , and the 
lefthanded quarks are placed in the (3, 1,2), then quarks can obtain masses from the following 
Yukawa couplings 

£ = - V3ggd m &Q Lj - V3g%u m e&*Q Lj + h.c. , (5) 

where $ and $' transform in the (3, 3, 2) and (3, 3, 2), respectively. 1 The fields $, $' can be 
decomposed into singlet and octet scalars under the unbroken SU(3) C . 

$(') = J_ H {>) + 5 a(/) T a (6) 

v3 

and A^' D = V3q^- D at tree level. Note that this model is different from the Manohar- 
Wise model since there are two distinct color-octet scalars, S a and S al . If the gluon-gluon 
fusion production of a single color-octet scalar proceeds through a top-quark loop, this 
mechanism will predominantly produce S al '. We do not wish to to study this model further, 
but mention it to show that the MFV constraint A^' D oc g^- D could emerge as a consequence 
of a symmetry of the underlying theory. 

The color-octet scalars can be produced via the pair production cross section, gg — > SS 



which proceeds through gauge couplings 16]. Constraints on pair production followed by 



1 Note that an additional scalar in the (1,1,2) is required to give masses to leptons. 
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decay to heavy quarks has been used to establish a lower bound on mj of about 200 GeV in 



Ref. [251 ] . Neutral color-octet scalars can also be produced singly via gluon-gluon fusion |2J 
This proceeds through loop diagrams containing quarks, of which the top quark gives by far 
the dominant contribution, and loops with color-octet scalars. The relative size of the top 
quark and scalar loop contributions is determined by rju and other parameters in the scalar 
potential. If these parameters are all taken to be of order unity then the top quark loop is 
the largest contribution. 2 In this case the production mechanism is very similar to that for 
a single Higgs, and hence threshold corrections are expected to be significant. At the LHC 
the production cross section for single S production is larger than pair production when the 
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mass of S is larger than 1 TeV 

When applying SCET to color-octet scalar production we first match full QCD onto SCET 
operators at the hard scale fih ~ ms, where mg is the mass of the color-octet field. SCET 
is formulated as an expansion in A ~ y /i s / fih, where /i s is the soft scale, ii s <y mg(l — z). 



The allowed SCET operators are constrained by SCET gauge symmetries 26j . At leading 



order in A, we find two dimension-5 operators with different color structures. Subleading 



operators will be constrained by the requirement of reparametrization invariance 27 
Our result for the factorized scattering cross section is 

f 1 dz - 

a(pp^SX) = rH(m s ,iif) / — S(m s {l - z), fi f )F(r/z, (i f ). (7) 

Jt z 

Here r = m 2 s /s where s is the center of mass energy squared at the LHC, H(ms,fi) and 
S(ms(l — z),n) are the hard and soft functions respectively, and F(x, /i) is the following 
convolution of the gluon PDF's: 

F(x,n f )= -^-f g/P (y 7 n f )f g/ p(x/y,n f ). (8) 

Renormalization group equations for the H(ms, /i/), S(ms(l — z),/j,f), and F(r/z,Hf) can 
be used to resum large threshold corrections as will be discussed below. 

In order to prove the factorization theorem in Eq. ([7]) we need to construct the SCET 
operators composed of the collinear gluons from the initial state hadrons, soft gluons, and a 
heavy color-octet scalar field. In the center of the mass frame, the incoming gluons from the 
two incoming protons are described as n and n-collinear fields where the light-cone vectors 
satisfy n 2 = n 2 = 0, n ■ n = 2. The lowest dimension operator with a single n-collinear gluon 
that is n-collinear gauge invariant is W\G^W n . Here G£" = G a ^ v T a is a SCET gluon field 
strength tensor, and W n is the collinear Wilson line 

W n (x) = P exp (ig J dsn- A a n {sn»)T^j . (9) 

Similarly, the lowest dimension n-collinear gauge invariant operator is W^G^Wn. Combin- 
ing W\G^W n and W^G^Wn into a Lorentz scalar and then expanding to lowest order in 



If the color-octet is a pseudoscalar only the top quark loop contributes. 
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A yields 

(wlG^W n ) (wiGn^Wn) = ~(wl[iD£,iD%W n ) (wi[iD nix ,iD mv )Wn) (10) 

V / a/3 \ / -yS g \ / a/3 \ / 70 

- -( e »")«,K) 7j + °w- 

Here a, (3, 7, and 5 are color indices in the fundamental representation and Bp is 



B^ = 

n 9 



1 n-VW^iD^Wn , (11) 



where the derivative operator P M returns the large label momentum and only acts on collinear 
fields within the brackets [•••]. It will be convenient to write the field Bp in terms of the 
Wilson line in the adjoint representation (i.e. with color generator (t a )b c = —if abc ). Defining 
Bp = B a ST a , B a S is given by 

Kf = ^g^W^Gl^ = mV^X", (12) 

where W n is the collinear Wilson line in the adjoint representation, and we used the relation 
= W^ a for the second equality. For the n-collinear fields, Bp and B^_ are identical to 
Bp and B°p_, respectively, after interchanging n and n. 

Finally, we need to include fields for the color-octet scalar. The strong interactions of 
this field are described by Eq. (j3J). At scales well below 777.5, the strong interactions of 
the heavy color-octet scalar simplify because the scalar is slowly moving. In the threshold 
region, the 5* a is produced nearly at rest (in the parton center-of-mass frame) and heavy 
particle effective theory techniques can be applied. We use a heavy scalar effective theory 
(HSET), similar to heavy quark effective theory (HQET). In HSET, the scalar momentum is 
decomposed into large and small parts: p$ = rrisv M + k 11 , where v M is the static four- velocity 
and represents fluctuations of 0(fi s ). In order for derivatives in HSET to bring factors 
of k^ rather than the total momentum, we use the standard rephasing trick to relate full 
theory and HSET fields : 

S a {x) = -jL= ( e- imsv - x S a Jx) + e imsv - x S* a {x)) . (13) 
V2m s V / 

The HSET Lagrangian is obtained by plugging Eq. (TL3]) into Eq. (j3J) and taking the large 
ms limit, in which the only surviving terms are those for which the phase factor cancels. 
We find 

/w = s* v a (v ■ iD s rs: - J-s: a (D 2 s rs c v , (u) 

2m s 

where is the four-velocity and the covariant derivative, D%, involves only soft gluons. The 
first term in Eq. (1141) gives the leading interactions, and the second term is suppressed by 
1/ms and so can be neglected. 

In our SCET-HSET operators, the soft gluons appear in the soft Wilson lines, 

y v (x) =Pexpf ig J ds v • A^sv^tA , (15) 
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where v M can be either n^, n M , or v^. These soft Wilson lines arise when we decouple the 
leading soft interactions from the collinear and heavy scalar fields by the field redefinitions 
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A a n » -> y^ b A^, w^ )ab -> y a n c w^ )cb (i6) 

71 ^ 71 71 ) ' ^71 ' v n 



n i 

ab ob o*a . -\-tba o*b 



After this field redefinition, the collinear fields and HSET fields do not interact with the 
soft particles. Note that after the field redefinition, £hset = S%*(v ■ id)S% + 0{l/ms), so 
the strong interactions vanish at leading order in 0(l/ms)- The heavy scalar 's interaction 
with soft gluons can be reproduced by a soft Wilson line. This replacement simplifies the 
derivation of the factorization theorem as we show below. 

Using B^, B^, S%, and the soft Wilson lines we can construct the effective Lagrangian 
for color-octet production at leading order in A, 

£scet = C s (fi)O s (fi) + C P (fi)Op(fi), (17) 

where the effective theory operators Os and Op are for color-octet scalars (Ss) and pseu- 
doscalars (Sp) respectively. Those operators have different color structures and are 



d abc 

y/ 2m s' 



Os = ^=(s* Sv yty(B^)\y n B^y, (is) 



• £abc 

Op = ^^4 v (s Pv ylr(B^yl)\y n B^) c , (is) 



where ej[ u = e^ pa n p n u /2. Note that the strong interaction Lagrangian for the pseudoscalar 
is the same as the scalar. Since the n-collinear, n-collinear, and the soft fields are decoupled, 
the renormalization of the both operators will be a simple product given by 



J S,P 



ZyZ n Zn, (20) 



where Z n and Z w are the renormalization factors for the collinear parts B^ and B^^, 
respectively, and Zy is the renormalization factor for the three soft Wilson lines. So the 
renormalizations of Os and Op are the same and do not depend on either color structure 
constants or the Lorentz structure. Therefore, below we only consider the renormalization 
of the scalar operators, however our results hold for pseudoscalar operators as well. 

Taking the matrix elements of C s{n)0 s{n) in Eq. (IT7|) . we find the scattering cross section 

a(pp ^ S S X) = ~J2[ d A q5{q 2 -m 2 s )5{P n + P^-q-p x ) 
s x J 

x\Cs(ms,fi)\ 2 (P n P^\Ol(ti)\S s X)(SsX\Os(fi)\PnP^, 

where P£ and P£ are the momenta of the incoming protons which are n-collinear and n- 
collinear, respectively. Because the final state X consists of n(rT)-collinear and soft states in 
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the partonic threshold region, it is possible to rewrite the final state summation as Ylx = 
Ylix n ^2x w ^2x s anc ^ the final state momentum p x = Px n + Px w + Px s • The momentum of 
the color-octet field is 



q = P n + Pn~ (p Xn + PX W + PX S ) =Pn+Pn-PX s 



(21) 



where p n = P„ — px n and Pn = Pn — Px- are momenta of the partons in the two incoming 
protons. Then the argument in the first delta function in Eq. ( |2~TT) is 

q 2 -m% = (Pn+Pn- PXsf ~ m \ ~ (Pn + Pn) 2 ~ ^PX S ' {Pn + Pn) ~ m| (22) 



s — 2i]s 1 ^ 2 — m 



s- 



where rj = p° Xs is the energy carried by final state soft particles in the partonic center-of-mass 
frame, and s = n ■ p n n ■ p^ = y\y 2 s, where yx j2 are the large (collinear) momentum fractions 
of the incoming partons, defined by yi — n ■ p n /n ■ P n and y 2 — n ■ Pn/n ■ P^. When the 
momentum q is on-shell, Eq. ( 1221) implies 



V 



s — rric. 
2s 



s(l — z) ms(l — z) ms(l — z) 
2 2zV2 2 



(23) 



where the last approximation is valid in the limit z = m 2 s /s = r/(yiy 2 ) — > 1. 
The cross section in Eq. ( I2ip can be factorized by first inserting the identity 



n ■ V 

dr]dy l dy 2 5{r] + id )6 ( y x - _ n 

n ■ P„ 



$ 2/2 



n ■ V 

n ■ PrT 



(24) 



where n ■ V is a label operator acting on n-collinear fields, n ■ V is a label operator acting on 
ra-collinear fields, and the partial derivative, id , acts only on soft fields. Then a(pp — > S$X) 
can be written as 



a(pp -> S S X) 



7T 



drjdyidy 2 5 ( rj - 



m s(l ~ z)\ \C 3 (m8,ii)\' 



2s 1 / 2 



(25) 



d abc d def 
X (P n Pn 

2m, 1 



{y v s Sv ) a {B^yi) b {y^ > 



nal 



SsX 



x(SsX 



s( v + tdo)(sM) d (B^]yiy{yj3^ [ yi ]) 



p P- 



The cross section a(pp — > SpX) is the same up to the color factor and the replacement 
C s (m s ,n) -> C P (m s ,fi). Then we use S%\S) = ^2l^e a and (S\S* b = V2^e* b to remove 
the color-octet scalar from the final state, where e a are color polarization vectors satisfying 
the relations, J2 pol e* a e b = 5 ab . Finally, we apply the completeness relation, 1 = \X)(X\. 
After these manipulations, we find that a(pp — > SsX) is given by 



a(pp -> S S X) 



7rd abc d def 



dyidy 



\Cs{ms,fj)\' 



P P- 

1 n 1 n 



x6 



m s (l - z) 



2sV 2 

^ )yl kd B^[y 2 }ytyl o B^ o [ yi 



P P- 
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To simplify the notation we have defined 

n-V 



S[Vi- 



n ■ P„ 



, B^ a [y 2 ]= 6[y 2 



n ■ V 
n ■ Ptt 



B 



-Lfj,,a 



The PDF for the n-collinear proton in terms of SCET fields is 



fg/P n ( x ) 



1 



dn ■ z 



-ixn-P n n-z/2 



2nxn ■ P n 
xn a nPg^(P n 



(27) 



(28) 



ah 



n ■ z 







1 



-n 



xn ■ P„ 
1 

x(n- Pr,) 2 



nPgT [Pn G a n ^W: b 5(xn -P n -n- V)W^ C G^ V 
n-V 



Pr, 



Pn 



B^Slr 



B ±,a 

nfJ. 



P 

- 1 n 



where P% is the momentum of the proton. Here, we have defined W^ 6 [z, 0] = W" !, (z)W| fc (0) 
and used Eq. (jl2jl for the third equality. A similar set of manipulations yields the PDF for 
the n-collinear proton, and since f g /p n (x) = f g /p w (x), we will drop the subscripts n and n in 
what follows. Therefore, after averaging over proton spins in Eq. ff26l) . we find 



(Pn\B^ l B^°[ yi ]\P n ) = 9T^ V ^^ fMvi) t 
The soft function S$ (Sp) for the scalar (pseudoscalar) production is defined to be 



(29) 



S s (m s (l - z)) 



3d abc d def 
40 



y^y^y^sfi - z + 2 l -^-)yi kd y^y^ 

V m.q/ 



£abc fdef 

S P (m s (l-z)) = J —^—(0 
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yfy^ m s( i-z + 2 l -^-)yl kd yl le y^ 



fmf 
n 

mf 





oV 



(30) 
(31) 



Here, the soft functions are normalized to 5(1 — z) at lowest order. 

Using these definitions we obtain the factorized scattering cross sections which are given 

by 

- 1 dz 



f dz - 

o{pp -> SiX) = THi(m s ,fJ.f) — Si{m s (l - z), n f )F(r/z, fi f ), i = S,P, 

J T Z 



where the hard coefficients Hi are 



H s {m s ,iif) 



57r|C , 5 (m 5 ,/i / )p 



H P (m s ,Hf) 



37r|Cp(m 5 ,/i/)p 



(32) 



(33) 



48 ' irjJ 16 

This factorization theorem is our main result. When we consider the renormalization group 
evolution effects, we will calculate Hs t p(ms, n) (Ss,p(iris(l — z), u)) at the scale \ih (p s ) and 
then evolve them to the factorization scale /if. At the leading order (LO) in a s , the cross 
section is 



a(pp - S,X) = rH^F(r) 



S,P, 



(34) 
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FIG. 1: One-loop renormalization of Os,p- Here the curly lines with the straight lines are n(n)- 
collinear gluons and the only curly lines denote the soft gluons coming from the soft Wilson lines. 
Double line denotes outgoing color-octet field. 

where H- ' are the hard coefficient at the lowest order and are equal to the scattering cross 
section at the Born level. 

Next we discuss the RGE evolution of the hard and soft parts and the resummation of the 
cross section. Once the evolution for the coefficient functions, Cs,p{n), is determined, the 
evolution equation for the soft function can be easily derived, since the evolution equations 
for fg/ p (x) are known. The evolution of the soft functions can be done in momentum space as 
in the analysis of Higgs production and Drell-Yan in Refs. 

119 

. We follow this approach 

in this paper. Alternatively, one can solve evolution equations for the moments of the soft 
functions and PDF's, and then take an inverse Mellin transform to obtain the resummed 
cross section. Resummed expressions for the moments of the cross section are given in the 
Appendix. 

To determine the one-loop anomalous dimensions of Os,p, we need to consider the Feyn- 
man diagrams in Fig. [1] as well as the wavefunction renormalization graphs. We regulate 
ultra-violet (UV) divergences using dimensional regularization and the infrared (IR) diver- 
gences by taking the external legs to be off-shell. It is then straightforward to extract the 
UV divergences and we find 



1 + 



o:. s 



N r 



47T£ UV L V£ UV 



2 +21,^ 



+ 



14 



mi 



+ m ) rif 



(35) 



where N c and rif are the number of colors and flavors, respectively. From Z$,p we obtain 
the anomalous dimension for Os and Op, 



ls,p 



7T 



NJ In 



mi 



7 7T 

H h i~ 

3 2 



3 



(36) 



Note that 



In Ar + i- = - In 



mi 



2 — m| — ie 



1 /i 2 



mi 



(37) 



so logarithms of both +m| and — m 2 s appear. This is because there are corrections coming 
from soft exchanges between the two initial state particles, similar to Drell-Yan, which give 
rise to logs of — m|, and soft exchanges between initial and final state particles, similar to 
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deep inelastic scatttering, which give rise to logs of +m|. From Eq. (I36I) we can infer the 
double logarithms in the 0{a s ) corrections to Ci(/i), i = S,P, 



(0) 



1 - 



47T 



log 2 



-m s — ie 



if 



+ 2 l0g2 



+ 



(38) 



where ... denotes terms without double logs. From this we see that if fi = ms, Ci(ms) 
gets a 7r 2 -enhanced contribution: Cj(ms) = C\ ■ (1 + C aols^s)^ 1^)- F° r the range of m 3 
considered in this paper, a s (ms) < 0.1. If « s = 0.1, this 7r 2 -enhanced correction increases the 
cross section by about 24%, and is half as big as the corresponding 7r 2 -enhanced contribution 
to Higgs production. Refs. [!, 0] argued that the 7r 2 -enhanced terms dominate the fixed- 
order corrections to Higgs production, and that these terms can be resummed to all orders 
by evolving the hard function from the scale m 2 H to the scale —m 2 H . They also showed that 
the leading terms exponentiate. In our case, setting /i 2 = — m| does not remove the factor 
of 7r 2 in the hard coefficient. The double logs vanish if 



fi 



3 7T(±l-i)/2 



III 



s- 



(39) 



but it is not clear that evolving to this complex scale will give a sensible resummation the 
7r 2 -enhanced contribution. Below we will calculate the i^-factor with the next-to-leading 
order (NLO) 7r 2 -enhanced correction. Even if the 7r 2 -enhanced corrections exponentiate, the 
NLO correction should be a good approximation to the resummed result since 1 + Caol s ^ '/4 
and exp(CUci!s7r/4) differ by less than 3% for a s = 0.1. 

The soft functions in Eqs. ( 1301) and (pTCl) can be computed perturbatively when /x s ~ 
777,5(1 — z) ^> Aq CD . The Feynman diagrams in Fig. [2] give us the 0(a s ) corrections to 
S s (m 3 (l -z),n) and S P {m s {l -z),n)\ 



S ( s\m s (l -z),n) 



7T 



1 



p (m s (l 



(40) 



is 



+ 



1 



-( 



uv 

— + l + 21n4- 

e uv mi 



E\jv v2 

1 



1 1 t* 2 

- + lnAr 



1-z) 



+ 4 



ln(l 



— + - In + - In 2 

4 2 m 2 s 2 m|j 



where the plus distributions are defined in the standard way. Note that there is no IR 
divergence in the sum of the real and virtual diagrams in Fig. |2] 28j. The IR finiteness of 
the soft function can be easily understood in SCET because the soft function is just the 
Wilson coefficient obtained at the second-step matching. 

To obtain the resummed cross section we employ the method of momentum space resum- 
mation developed in Refs. l29l . |3Cj • The resummed result can be written as 

f 1 dz 

a(pp^SiX)=r —Vi(z,m s ,fif)F(T/z,Hf), i = S,P 



where the resummation function Vi(z,ms, fif) is given by 

7. 

Vs t p(z,m s ,iif) = H s ,p(ms, Hh)U(fi h , fi s , fi f ) 



1 - zf-^ 



(41) 



(42) 



10 



(a) 



(b) 






FIG. 2: One loop corrections to the soft function. The dashed line represents the cut. The diagram 
(a) and its Hermitian conjugate (b) describe the virtual soft gluon radiation and the diagram (c) 
denotes real soft gluon radiation. 



Here Ss^pid^, is defined in terms of the Laplace transform of the soft functions, and 
the evolution function U(fih, fi s , /^/) is a product of terms obtained from evolving the hard 
function to the scale Hh and the Laplace transform soft function to the scale fi s . To NLL 
accuracy, we find 



In U(nh, ii s , fif) = In 



^SU NLL {m s , Us) + -5- In — -, r 

Po oi s {m s ) 



a. 



where S77nll is 



-S^nll (/•*!, ^2) 



47T 



a a (//i) 



1 In r 

r 



+ U"^ (1 " r+lnr)+ ^ ln 



+ 7T ln , 



(43) 



(44) 



where r = a s (/^2)/a s (/xi). The parameters Ai, A2, Bf , and B\ appear in the anomalous 
dimensions of the hard and soft functions, and are given in the Appendix. The parameter 77 



is defined in terms of an integral over the cusp anomalous dimension (see Ref. 29|) and in 
our calculation 77 = (Ai/f3 ) ln(a s (/x/)/a s (yU s )). In our case, r\ < since fi s < fif and hence 
the integral in Eq. (l4"Tj) is singular. The integral is then defined by analytic continuation 
from positive 77. We will choose [ih = = ms- F° r this choice of \ih there are no large logs 
of rrig/fi 2 in H St p{m s , fih) and H Sj p{m s , fih) = Hf p (m s , fit) to NLL accuracy. In order to 
resum logarithms of 1 — z we should choose the scale fi s = ms(l — z), however, this will 
lead to divergences in the z integral as the running coupling will cross the Landau pole as 
z — > 1. Practically, it is simpler to choose /x s to be a scale parametrically smaller than \i h . 

We first present our results in terms of a i^-factor, defined as the ratio of leading order 
and NLL resummed cross sections, which is given by 



/ 1 dz I 



1 dz 



(45) 

where i^o is a convolution of PDFs at LO. This result is universal in that it is independent 
of the NP model. 0(as{ms)) corrections to the hard coefficient can depend on the NP model 
but this beyond the accuracy we are working. For our numerical calcuations, we use the LO 
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V7 =14TeV 
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5(1(1 



1000 
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M s (GeV) 
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3000 



FIG. 3: K-factor for the single color-octet production at LHC where yfs = 14 TeV. The straight 
(dashed) line denotes NLL evaluation with (without) 7r 2 -evolution. 



a s , setting a s (M z ) = 0.1205 and m t = 170.9 GeV. For the gluon PDF's we use CTEQ5 at 
NLO 31]. In order to determine /z s , we follow the procedure of Ref. [2J and calculate the 
convolution of the one-loop expression for Si(m s (l — z),fi) in Eq. (T4"0l with F(t/z,[i) in 
Eq. ([7]). The scale fi s is chosen so that the higher order corrections to the soft function are 
under perturbative control. This is accomplished by specifying //f and fi 1 /: The scale /4 is 
defined by starting with /i s = /i^ and lowering /i s until the 0(as) correction is less than 15%. 
The scale [i 1 / is chosen so that the one-loop correction is minimized. We use the average 
(/if + fi 1 /) /2 in Fig. [3j The solid line is the result for the i^-factor with the 7r 2 -enhanced 
correction included. The .ff -factor varies from about 2.4 for ms = 500 GeV to about 3.6 for 
ms = 3 TeV. As expected, the resummation of threshold corrections significantly enhances 
the cross section and becomes more important as ms increases. The dashed line in Fig. [3] 
is the result without the 7r 2 -enhanced correction. This correction increases the .fT-factor by 
25% and is independent of mj. In Fig. HJ we show the variation in the prediction as fi s is 
varied between fi ! s and fi ! s ! . The uncertainty from the choice of fi s is ±15% for ms = 500 
GeV, and decreases with increasing ms- The variation with the choice of /x/ is also shown in 
Fig. HI The sensitivity to the choice of fif is greater and introduces an uncertainty of ±25%. 
The dependence on the scales fi s and /i/ should decrease when higher order corrections are 
included. 

In Fig. [51 we show our calculation of the color-octet scalar production cross section in the 



Manohar-Wise Model 16]. In this model, the two real components of the complex color- 
octet scalars are denoted and S®, where is a scalar and S® is a pseudoscalar if T]a is 
chosen to be real. The LO calculation of their production cross sections from Ref. [24J are 



the dashed lines in Fig. HJ and our NLL results are the solid lines. At ms = 1 TeV, we obtain 
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FIG. 4: Scale dependences of the K- factor. 



yfs =14TeV Vs~=14TeV 




500 1000 1500 2000 2500 3000 500 1000 1500 2000 2500 3000 



M s (GeV) M s (GeV) 

FIG. 5: The scattering cross sections employing Manohar-Wise Model at the LHC. In the both 
plots, the straight (dashed) lines denote the results at NLL (LO). 



0"nll(pp - > SrX) = 57 fb and — * SjX) = 73 fb. We have fixed the parameters 
T)u — 1 and A^s = 1 as in Ref. [24]. The NLL results are almost 3 times as large as the LO 
results, which are cr^o(pp — ► SrX) = 21 fb and <Tlo(pP — > SjX) = 26 fb. 

In summary we have used SCET to derive a factorization theorem for color-octet scalar 
production at the LHC. The factorization theorem can be used to resum large threshold 
corrections which have a significant impact on the total cross section. It is universal in 
the sense that all details dependent on NP models are encoded in the Wilson coefficients. 
The factorization theorem is similar to Higgs production, however, some details are different 
because the final state particle is colored. Because there are both soft exchanges between 
initial state partons as well as between partons in the initial and final states, the structure 
of double logarithms and corresponding 7r 2 -enhanced corrections is different. We obtained a 
resummed calculation of a(pp — > SX) to NLL accuracy. The resummed cross sections are 
2-4 times larger than the LO cross section, depending on the mass of the color-octet scalar. 
Uncertainties from varying [i$ and \Xf in these calculations are ±15% and ±25%, respectively. 
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NNLL log resummation and higher order perturbative corrections will be required to reduce 
scale dependence of the resummed cross section. Further development of the factorization 
theorem to account for scales besides m s and 777,5(1 — z) maybe required if color-octet scalars 
are actually discovered. For example, precision measurements of the mass may require taking 
into account the width of the color-octet scalar, as is required for determining the top quark 



mass 



321. 
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APPENDIX A: LARGE N RESUMMATION IN MOMENT SPACE 

Here we present our results for the resummed cross section in moment space. All the 
results below are taken in the large N limit. For the renormalized soft function Sn we find 

S£p(aO = C z N - x S s Am S {l-z)^) = l + ^Ar c (iln 2 ^ + iln^ + l + ^)+0(a s 2 ), 

(A1) 

where N = Ne 111 . From Eq. flAlj) , we notice that the choice \x = rris/N minimizes the large 
logarithms. The //-independence of cross section implies the following RGE: 

0") = (27? - 2Re[ 7< ])TO, i = S,P, (A2) 

where 7^ is the well-known Altarell-Parisi evolution kernel for the gluon PDF in the moment 



space, 



^ = -c 



71 



' A [2\nN-(^-^)\+0(al). (A3) 

From our results, Eqs. (j3SJ, flAip . and (|A3[) we can easily see that Eq. (IA2I) is satisfied to 
first order in a s . If we take the moments of a(pp — > SiX) in Eq. (!32l) . the result is 

a N (pp - SiX) = j\rr N - l cj{pp - S t X) = ^(m s , VsW(» s , ///)[/^(///)] 2 + O (J-) , 

(A4) 

where we identified \lk = ms- Here we employed the two-step matching: the hard coefficient 
Hs,p( m s) & t the scale m s is evolved down to the soft scale /i s and then the soft function 
S(fjL s ) obtained at fi s can be evolved to the factorization scale \if. This is equivalent to the 
scaling evolution realized in Eq. (142 j) . where the hard and soft function are evolved from ^ 
and fi s to jJLf respectively, but the soft function's renormalization behavior compensates the 
evolution of the hard function from fi s to fif. So the exponentiated matching coefficients 
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Cs,p and Sg P are given by 



13 d^i . ■ 
— 1s,pW 



Slp(»s,»f) = S^s)e- 12 ^ = S»(ii s ) exp 



-2 



A'/ 



7» 



Therefore Eq. (1A4I) can be rewritten as 

= fli(m s ) exp [eK/i/jJ^y 2 , 



(A5) 
(A6) 



(A7) 
(A8) 



where we set fi s = m s /N in the second equality, and then the exponential factor Q s can be 
expanded as 



Q(ms,fif) = kiNgP + 9s\ms,Hf) + a s (m s )g ( s' \ms,Hf) + 



(2), 



(A9) 



Here each of the coefficients g$\ i = 0, 1, 2 are correspond to the resummed results at LL, 
NLL, and NNLL accuracies, respectively. 

For the computation of the exponentiation factor up to NLL accuracy, we need to expand 
Re[75 ! p] and 7^ up to second order in a s 



Re[ 75 ] = Re[ 7P ] = [a, In ^ + 2?f ] - (g) A 2 In ^ + 0(a 2 s ), (A10) 
AilnJ^-Sf] + (^fA 2 \nN 2 + 0{a 2 s ), (All) 



An 



where a s hi(fi/ms) ~ a s In A" are treated as 0(1), and the coefficients of the large logarithms, 
Ai denote the coefficients of the cusp anomalous dimension 33j. In the above equations, Ai, 
B£, and Bf were already given in Eqs. ® and (|A"3| . and A 2 is 8A" c [(67/18-7r 2 /6)A' c -5r2 / /9] 



33 



After a brief calculation using Eqs. (1A10I) and (lAllj) . we find 



(0) 



9f 



1 

A) 
, PxAi 



2A + (1-2A) ln(l-2A) 



(Sf - 5f ) ln(l - 2A) + 2A X X In 



2 



PI 



(A12) 

2A + ln(l-2A)l (A13) 



4A + (2 + ln(l-2A))ln(l-2A) 



where A = a s f3ohiN/(A7r), and /?o,i are the first two coefficients of the QCD (3 function. 
Here note that we have set [if ~ ms. However, if we choose /i/ as fif < fi s ~ ms/N, the 
logarithm hi(fj,f/ms) should be power-counted as 0(lniV). 
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